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A dilute monopole gas explains, in quarkless gluodynamics, the small ratio 5 be- 
tween the square of magnetic screening mass niM and spatial Wilson loop tension. 
This ratio is 0.0895 for T = to 0.0594 at T = oo for any number of colours with 
order N~'^ corrections and equals up to a numerical factor of 0(1) the diluteness. 
The monopoles have a size (jv/ = ™jv/ ■ The GNO classification tells us they are in 
a representation of the magnetic SU(N) group. Choosing the adjoint for the dilute 
gas predicts the k-tensions to scale as k{N — k) , within a percent for high T, and 
a few percent for low T for the seven ratios determined by lattice simulation. The 
transition is determined by the transition of the dilute Bose gas at Tc = O.lYAmM , 
and the transition is that of a non- or nearly-relativistic Bose gas. 

1 Introduction 

An ancient idea in QCD ^ is that monopoles are responsible for flux tube formation through 
a dual superconductor mechanism. Another mechanism, that of the "Copenhagen vacuum" of 
the early eighties ^, proposes macroscopic Z(N) Dirac strings or Z(N) vortices. 

In this talk I will discuss a specific model of the first type"*, that works surprisingly well 
in its most direct applications. A straightforward way to see its workings is to start from the 
plasma phase. That sounds at first sight self-defeating as there is no confinement in this phase. 

But it will turn out that precisely the absence of confinement renders the detection of these 
monopoles, or perhaps more appropriately, magnetic quasi-particles, so straightforward. We 
know that electric quasi-particles, the gluons, are approximately free at very high temperature. 
This follows from the Stephan-Boltzmann law, which is verified up to 15% on the lattice, at 
T ~ ATc- Note, though, that the interactions at this temperature are still quite strong. 

We may guess that the same happens to the magnetic monopoles that were condensed in 
the cold phase. 

What kind of monopoles can we expect? In the absence of any spontaneous breaking the 
GNO classification'^ applies. For a theory with only Z(N) neutral fields, like the adjoint, the 
magnetic group is SU{N). So we can have monopoles in the fundamental, adjoint or any 
representation we like " . 

One of the consequences of the idea of a monopole condensate is the screening of the 
magnetic Coulomb force between two static magnetic sources. This is a straightforward gener- 
alization of the screening of the Coulomb force between two heavy quarks. There is now ample 
affirmation of screening * from lattice simulations. The magnetic screening mass rriM starts 
out at r = and equals there the lowest glueball mass. It stays constant till about 2Tc, 
and then starts to scale with the temperature like g^T. 

It is reasonable to think of the monopole as having a size of about the screening length 
Im — iTT-^i ■ We will assume that its size is much smaller than the inter-monopole distance. 
That, together with the choice of multiplet, will fix the ratios of spatial string tensions. The 
reader who is only interested in how one computes these ratio's should read sections 2,5 and 6. 

"In the GNO analysis the range of the monopoles is infinite. That is why quarks imply a smaller magnetic 
group, SU{N)/Z{N), because of the Dirac consistency relation between electric and magnetic charges. Also the 
existence of the monopoles was not proven. See e.g. ''' . 



The remaining sections give a qnalitative idea of how flux loop averages dep(^n(i on the 
representation. The main conclusion is that the loops depend only on the N-allity. Moreover, 
in our model only the fully anti-symmetric representations are a practical tool to determine 
numerically the tension. 

2 Electric flux loops 

In this section we discuss the behaviour of electric flux loops in the plasma, as it will explain 
the basic features that are permeating the talk. 

2.1 QED 

Consider a plasma of ions and electrons. We will take the ions to have the same but opposite 

charge of the electron. 

Suppose we want to compute the electric flux going through some large (with respect to 
the atomic size) closed loop L with area A{L). Normalize the flux = Jg dS ■ E by the electron 
charge e and define : 

F(L) =expi27r^/e. (1) 

Of course, at T below the ionization temperature no flux would be detected by the loop, 
because there are only neutral atoms moving through the loop. 

Let us now raise the temperature above Tionisation- What will happen? Both electrons and 
ions are screened. For simplicity we will take the ions to have the opposite of one electron 
charge. 

We are going to make the following simplification. The charged particles are supposed to 
shine their flux through the loop if they are within distance Id from the minimal area of the 
loop. This defines a slab of thickness 21d- 

Of course if we plot |^/e| as function of the distance of the particle to the loop you find an 
exponential curve with the maximum 1/2 at zero distance. For the sake of the argument we will 
replace that curve by a theta function of height 1/2 and width 2/d- If one wants to do better 
one has to deal with infinitesimally thin slabs, and integrate over the thickness. This yields an 
effective thickness of 1.64282... Here we will keep the factor 2. The result is parametrically 
the same as the one we will derive keeping the simple minded method. 

Then one electron (ion) on the down side of the loop will contribute +l/2(— 1/2) to the flux, 
and with opposite sign if on the up side of the loop. That is: V{L)\onc charge = ^1- This result 
is independent of the sign of the charge! The plasma is overall neutral, the loop is sensitive 
only to charge fluctuations. For I charges inside the slab the flux adds linearly and we find: 

v{m = {-iy. (2) 

Assuming that all charges move independently, the average of the flux loop V{L) is deter- 
mined by the probability P{1) that / electrons (ions) are present in the slab of thickness 2iD 
around the area spanned by the loop. Taking for P{1) the Poisson distribution jf(r)'exp— F- 
l is the average number of electrons (ions) in the slab- we find for the thermal average of the 
loop due to both ions and electrons: 

(f(l))t = Y,PimL)\i = Y^pm-y = exp-4r. (3) 

; I 
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Now / = A{L)2lr)n{T) , so the electric flux loop obeys an area law ex.p—p{T)A{L), with a 
tension 

piT) = SluniT). (4) 

This area law distinguishes the behaviour of the loop in the plasma from that in the normal 
unionized state. 

The flux loop has a very useful alternative formulation. Introduce the charge density jo- 
Then you can write flux loop as a gauge transformation exp ( — f / dx{E.V -\- jo)u)i)- Wi a 
gauge function that falls off fast atspatial infinity. But it has a discontinuity 27r, when crossing 
the surface. So the border of the surface is of course just a closed Dirac string. 

If LOi has no discontinuity, integration by parts will give us the Gauss operator — V.S + jo 
in the exponent. But the discontinuity generates the extra surface term Because the flux 

$ is gauge invariant the two terms commute. So the gauge transform factorizes in a factor 
containing the flux and a factor containing only the Gauss operator. The latter becomes unity 
on the physical subspace and only the flux operator stays 



2.2 Electric flux in gluodynamics 

Is there an SU{N) generalization of the QED case? In fact yes. The closed Dirac string in QED 
is replaced by a Z(N) vortex of strength fc, the 't Hooft loop?. We introduce in the Lie algebra 
of traceless hermitean N x N matrices a basis for the Cartan subalgebra, the {N — 1) x {N — 1) 
dimensional subspace of diagonal matrices. This basis of — 1 diagonal matrices Yfc is chosen 
such that in exponentiated form exp ( — i2'jTYh) it gives the N—1 centergroup elements expifc^. 
A simple choice is 

Yk = ^diag{N - k, ,N-k,-k, ,-k). (5) 

The entry N — k comes in k times, and —k comes in N — k times, so the trace is 0. 
The flux operator becomes in this notation: 



Vk{L)=expi^ j dS.TrEYk. 



(6) 



It does correspond to the vortex operator of 't Hooft with strength k 
The next question is: does the gas of deconfined gluons induce an area law in this operator? 
The answer is yes, and the reasoning is as before. The charge of a gluon with respect to the 
charge Yk is found from the form Yk takes in the adjoint representation. This is easy: we 
have charge ±1 like in QED, but unlike in QED we have 2k{N — k) gluon species in the gluon 
multiplet with such a charge. All other gluons have vanishing charge so do not contribute to the 
flux. Since we take the gluons to be statistically independent, the charged ones all contribute 
a factor to (Vfc), and this factor is the same by SU{N) symmetry: exp—ilon. 
Conclusion: the expectation value of the loop is 

{Vk)=exp{-pk{T)A{L)) (7) 

with the tension 

Pk=AlDnk{N -k). (8) 



''Normalization is such that [B? (0) , (5)] = 5°-'''5j^k5{x), E = E°-Xa,{\a,\b\ = ifa,b,Ac,fa,b,cfd,b,c 
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So this is the k-scaling law for the electric fliix loop. It does obey large N factorization: 



Pk = kpi. 



(9) 



But it has 1/A'' corrections and not, as perhaps expected in gluodynamics, 
This computation is corroborated by low order perturbation theory for the loop. The 
expectation value of the loop can be computed from a tunneling process between adjacent 
vacua of an effective potential with a Z{N) symmetry This potential has, because of the 
Z{N) symmetry, degenerate minima in the centergroup elements. Only at cubic order in the 
coupling g the k-scaling starts to deviate slightly Simulations in S'C/(4) show k-scaling up to 



The effective potential is computed in a background that serves as profile orthogonal to the 
area. The minimal effective action of this profile gives the tension pk, thereby tunneling from 
say 1 to expz27ryfc. The profile is present in the propagators and renders the usual double line 
representation of the gluon propagators invalid, hence the 1/N corrections. The typical width 
of the profile around the area is the Debye screening length Id- 

2.3 Junctions and highest weights 

It may have struck the reader that we used Y/j, instead of kYi in the flux formula cq.(6). The 
reason is that the tunneling process from to Yf^ passes only one mountain. On the other hand 
tunneling from to fcYi passes through k mountains, each of which interposes between two 
subsequent Z{N) vacua and contributes pi. So at first sight we find kpi from this operator. 
However, there must be a lower state pk, as the following argument shows. Formally the 
operator 



is obtained by coalescing the product of k operators V\{L). There are two relevant stages here: 
one where all the loops are still more than the screening distance away from each other. Then 
every loop will form its individual profile and wall, and give a factor exp —piA{L). The second 
stage is when all loops are within screening distance. Then the individual walls emanating 
from each loop will form a junction at a typical distance Id- At this junction k unit walls 
merge with the k-wall, and the k-wall views the ensemble of the k unit discontinuities as one 
discontinuity of strength k, on the scale of the screening length. The junction coupling involved 
here is non-perturbative. It may be small, certainly in our model of approximately free gluons. 
Hence the true pj. may be difficult to disentangle numerically 

In ref. the case k = N in eq.(lO) is discussed. The operator Vk=N creates a state with 
energy Npi and involves the special junction where k = N. This junction has only N incoming 
k=l walls, and no outgoing N wall. That explains the periodicity A'^ of the tension pk- 

In addition we can imagine junctions where two of the incoming unit walls are replaced by 
one k=2 wall. That junction corresponds to the matrix (fc — 2)Yi + Y2. 

A general jimction in SU(N) will be specified by a set of non-negative numbers wi,W2, , wn-i- 

So the constraint on those numbers is: 



In our two previous examples we had wi = k all others zero, and wi = k — 1,W2 = 1 , 
all others zero. The corresponding matrix is then Xl;=i wiYi and the reader will recognize 



2Te 




(10) 



JV-1 




(11) 
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Figure 1; A Young tableau with uii = 0, 102 = 1 = W3, 104 = = U15, u>6 = 1, Hn = 5^2 + ^3 + 5^6 and k = 10. 

this matrix as being the highest weight Hji of the Young tableau characterized (see fig.l) by 

the numbers wi, wn-i^- These Young tableaux refer to the representations of the magnetic 

GNO group. 

Also this general junction J{k; {w}) is expected to be small for the same reasons. Especially, 
because of factorization, eq. (9), it vanishes in the large N limit. The conclusion is that only 
for the partitioning Wk = 1 , all other w's zero, wo obtain a direct coupling to the k-tcnsion, 
and no junction is needed. All other charges will give states of non-interacting k-tensions. 

Something very similar happens for the various representations involved in the Wilson loops 
and we will come back to this in detail in section 3.3. 

3 Some basic facts in quarkless Yang-Mills 

In this section wc gather a few facts, partly empirical, from lattice simulations, and partly from 
simple physics arguments. 

3.1 Magnetic screening 

A feature that sets gluodynamics apart from QED is the appearance of magnetic screening. 
When we put two heavy monopoles far apart at a distance r we find a Yukawa type potential: 

V{r) = - exp -niMr. (12) 

A simple argument shows * that the magnetic mass is the mass of a scalar excitation of a 
Hamiltonian in a world with two large space dimensions and one periodic mod 1/T. The space 
symmetries in such a world are SO {2) x P x C x R. We have the rotation group in the two large 
dimensions, C is charge conjugation, P is parity in the large 2d space, and R is the parity in 
the periodic direction and a state is written as J^*^^. The state corresponding to the magnetic 
mass is Oljl^. As T — > the rotation group becomes SO{3), and R and P become related 
through a rotation, so at T = the 0+"'" glueball mass results. The Debye mass, screening the 
electric Coulomb force between two heavy quarks, is related to the 0't~ state, and disappears 
below Tc- Below Tc Z{N) symmetry is restored and the string tension appears, as a state 
transforming non-trivially under Z{N). 

In contrast, the magnetic mass does not feel the Z{N) transition. At temperatures where 
T dominates all scales, the magnetic mass starts to scale like g^T, as schematically shown in 
fig. 2. 
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3.2 Spatial tension 

The spatial tension is defined througli a spatial loop L, on which an ordered product V cxp i § ds.An 
lives. R labels the representation for the vector potential we have chosen. The spatial Wilson 
loop is then defined as the trace over colour: 

Wfl(L) = TrV exp i j> (Is.Ar. (13) 

The path integral average gives an area law: 

{Wii{L)) = Ce^^{-(jR{T)A{L)) (14) 

for all temperatures. At zero temperature the spatial tension equals the string tension ° from a 
time-space loop because of Euclidean rotation invariance. But as function of temperature they 
behave differently. The string tension suffers a correction at low temperature of — f due to 
the excitations contributing the Luescher term. Above Tc it vanishes. The spatial tension stays 
flat, see fig. (2). Like the magnetic mass it does not feel the Z{N) transition and starts to scale 
like [g^TY well above the critical temperature. The dimensionless ratio of the spatial tension 
in the fundamental representation and the square of the magnetic mass is small, on the order 
of a few percent. In two extreme cases the ratio 5{T) is accurately known from simulations : 
at T = cxD * in the 3d theory, and at T = In between there; arc simulations that 
are consistent with a slowly varying 5{T). It is very important to have this verified with more 
precision. 



3.3 Dependence of the tension on the representation, and N-allity 

There is widespread agreement on the dependence of the tension and the string tension on the 
representation R. It is only its N-allity k that counts. If we think of R being built up by / 
fundamental and / anti-fundamental representations, the N-allity is just the difference: 

k = f-f. (15) 

To get a more precise idea, let us imagine the Wilson loop is formulated in a periodic box, 
with the 3 space dimensions being of macroscopic length L. The fourth dimension is of length 
1/T. The loop is now replaced by two Polyakov lines of opposite orientation in the x-dircction. 
They are a distance r apart- and carry the reprsentation R. The expectation value of the loop 
is then parametrized as: 

{Wr{L)) = Ce^v-LVrir). (16) 
If the distance r is very small, asymptotic freedom tells us: 

^ ^ Airr ^ ' 

where C2{R) is the quadratic Casimir operator. 

In three dimensions the Coulomb force is logarithmic. 

The Casimir operator can be related to our NxN matrices Yfe if R corresponds to a Young 

tableau, with the first row having wi more boxes than the second. The second row has W2 more 
boxes than the third, and so on. By definition, the numbers Wk are never negative. 

°In what follows we reserve the name string tension for the time-space loop. 
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Figure 2: Magnetic mass niM and tension a as function of temperature, schematically. mjv/(0) = mQ++ and 
from ref. (13): Tc = 0.174 mQ++. The temperature Tq < mQ++ is where the de Broglie thermal wave length 
becomes equal to the magnetic screening length. For the calculation of the tension it is below Tq that quantum 
statistics applies, above classical statistics applies as in section 5. 



There is a one-to-one correspondence between Young diagrams and irreducible representa- 
tions of SU(iV). 

Define the highest weight of R as 

Hn^J^WkYk (18) 

k 

as we aheady alluded to at the end of section 2.2 . An example is shown in fig.l. 
Then, if y = Yk : 

C2 (R) = i {TrHl + 2TrYHR) . (19) 

If the distance is long enough ( a somewhat ambiguous criterion!) the string regime sets in 
and we expect: 

V{r)=ak{T)r. (20) 

The tension only depends on the N-allity k. String formation between the two sources 
renders this dependence very plausible. The tension ai of the string formed in between two 
fundamental sources can form a bound state in the string with tension (72. We can go on like 
this and the question is, what is the dependence on fc? It should obey certain a priori criteria. 





® 











Figure 3: Young tableaux with k=2 

For large number of colours we expect factorization 

ak = kui. (21) 
Prom SUSY experience one also expects in that same limit: 

ak = kai+0{l/N^). (22) 
This is discussed in depth by Misha Shifman in this volume. 

3.4. Flux representation of the Wilson loop 

If wc want to proceed with the Wilson loops, as we did with the 't Hooft loops, we need clearly 
a representation similar to eq.(6). 

Remember we discussed in section 2.3 a general formula for the electric flux loop with Y}. 
in eq.(6) replaced by the highest weight but still with N-allity fc:VRexp (^ / dS.TrEHu). 

Of course a natural guess is to take this electric flux formula, replace E by B, a = g^/4:'jT 
by 1/a to get: 

Wr " =" exp {ig j dS.TrBHR). (23) 

This formula cannot be true^^. But it is true in a dynamical sense: consider the limit of a 
path integral average in the 3d theory with an adjoint Higgs, the electrostatic QCD Lagrangian. 

For a given representation R with highest weight Hr with stability group S (i.e. SHrS~^ = 
Hr) one flnds: 

{Wr) =< j Dfiexp {ig j dS.{B„ - -^f^'^Drn, A Dn,)Tr^nHRn'^) > . (24) 

with Ha = Tr^^HRfl^ the Higgs field's angular part that parametrizes the coset SU{N)/ S{Hr) 
The second term in the exponent is the source term for the monopoles in the sense that it carries 
no long range effects in the original Higgs phase. The average of the r.h.s is the 3d Yang-Mills 
average. The integration is over all gauge transforms with an arbitrary number of hedge hog 
configurations. The r.h.s. is obtained by taking the average of the magnetic charge operator 
in the 3d Higgs phase characterized by Hr. Then the the VEV of the Higgs field is let to zero, 
which introduces the fluctuations over n^, the angular components of the Higgs field. Finally 
one lets the mass of the Higgs become infinite, which siippresses the radial integrations. 

On the other hand the integrand of the r.h.s., without the brackets, was shown by 
Diakonov and Petrov to equal Wr, by one dimensional quantum mechanics methods. It involves 
a limiting procedure which is procured in the path integral by starting out from the Higgs phase 
and moving to the symmetic phase as described above. With this specification we will use eq. 
(23) without quotation marks in what follows. 

So we need a stable Higgs phase for the limiting procedure. The question is then: what 
are the stable Higgs phases? The answer is from simulations: only those that have a stability 
group S = SU{1) X SU{N - I) x U{1)'^^. 
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Summarizing: it is possible to have a magnetic flux representation in a path integral average 
for Wilsonloops Wr with special highest weight Hr. The highest weight should commute with 
a subgroup of the form SU{1) x SU{N—l) xU{l). These invariance groups define the only stable 
Higgs phases in 3d. The monopoles in these phases are partly screened in the strongly coupled 
sectors ( the non-abelian factors in the stability group). The SU(5) case is nicely described by 
Coleman in ref.^. After the limiting procedure they are the fully screened monopoles in the 
symmetric phase. 

To get a feeling, look at the k=2 representations of SU{N) in fig. 3). The highest weight 
of the totally anti-symmetric representation is computed from eq.(5) and eq.(18). It is Y2, 
and SU (2) x SU {N — 2) x U{1) is its invariance group. So it has a flux representation Wr = 

expig J dS.TrBY2. 

The symmetric k=2 representation has highest weight 2Yi. It has invariance group SU (N — 
1) X U{1). So its flux representation is Wr = expig J dS.TrB2Yi. 

4 Monopoles cind the magnetic group 

In this section we briefly touch on a medley of problems related to magnetic screening. 

The generators of the magnetic group as advocated by Goddard et al. ^ are not known in 
terms of the Yang-Mills potentials. In our opinion the monopoles are collective excitations of 
the magnetic gluons and can be much better called magnetic quasi-particles. They are bound 
states of magnetic gluons. Their size is supposed to be Im, the magnetic screening length. 

The question is then: what representation of the magnetic group is realized by Nature 
for these dynamical quasi-particles? We have tried in the next section the simplest ones, 
the fundamental and the adjoint. The adjoint has the advantage that it is compatible with 
quarks. Compatible means the Dirac consistency condition is fulfilled. As already mentioned 
in the introduction, the consistency condition is stricly enforced in the absence of screening. If 
screening is present a famous example due to 't Hooft ^ tells us that screening renders the Dirac 
condition less restrictive than naively thought. At any rate the adjoint is clearly favoured by 
simulations of the Wilson loops, see section 6. 

Amusing, though perhaps academic, is the observation that the flux representation of the 
't Hooft loops introduced highest weights Hr. These should correspond to representations 
of the magnetic group. If we could implement the loop with a given weight R by means of 
a magnetic gauge potential we would have an alternative expression for the 't Hooft loop. 
Once the magnetic group becomes a gauge group with gauge potentials one has an alternative 
for the Wilson loops as well, to wit, as Z{N) discontinuous gauge transformations like the 't 
Hooft loops. The magnetic gauge group provides long range excitations in the low T phase 
in contradiction with the full magnetic screening. They should be obliterated by a Higgs 
mechanism that breaks fully the symmetry. And for that our adjoint multiplet is not a good 
candidate. It will always leave some some subgroiip unbroken as we discussed in the previous 
section. A fundamental multiplet would fully screen, but is excluded by the lattice data. 

5 Predictions for the Wilson loops 

Once we are given a dilute gas of screened monopoles, we only need to specify the representation 
of the magnetic group for our monopoles. Then the calculation of the tension is done with our 
flux representation for the average of the loop, eq.(24). But in practice we can use the simple 
formula eq.(23). The reason is that we put the contribution of the monopole gas in by hand. 
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That takes care of the singular gauge transforms in cq.(24). The remaining regular ones can 
be dropped because we compute something gauge invariant. 

So the average of a k-loop in the totally anti-symmetric representation is computed from 

{Wk) = (expi5 j dS.TrBYk). (25) 

For the dilute gas in the adjoint representation the computation is not different from the 
one with gluons in section 2.2. The adjoint monopoles have a magnetic charge equal to 
or 0. Only the former contribute and their individual contribution to the loop in eq. (25) is 
— 1. The diluteness, and classical statistics (as the thermal wave length is and the typical 
interparticle distance is {g'^T)~^ this is justified, see also fig. 2.) then give for every charged 
species the same contribution exp — 2/m'^Mj being the common density of a given species. 
As there are 2k{N — k) charged species in the adjoint, the total results in the k-tension being: 

(Tfc = Ak{N - k)lMnM- (26) 

One can do a similar calculation for monopoles in the fundamental representation. The 
counting goes now as follows. Recall eq.(5), = j^diag{N — k, ....N — k, —k, .... — k) with trace 
zero. The Yfc charge of the highest k components of the N components of the column spinor is 
given by '^'^^^^ ■ The remaining N ^ k components have charge —j^- 

It then follows from the use of the Poisson distribution that the flux of a given component 
is contributing cos(7r^^^^) or cos(7r-^). Taking into account the degeneracies the final result 
for the tension becomes: 

cTfe = 4lMnM(iV - k cos(7r ^ ^ ^ - {N - k) cos{n—)). (27) 

Both tensions give ratios ak / fxi that behave like k for large N and finite k. This factorization 
is expected for a k-loop tension caused by screened particles'' . But for large N and k = N/2 
the ratio for the fundamental multiplet is a factor 2 larger. 



6 Comparison to lattice simulations 

We have been discussing a model at very high temperature. Hence it is tested in 3d lattice 
simulations. The ratios found for the totally antisymmetric irreps are close — within a 
percent for the central value — as far as the adjoint multiplet of magnetic quasi-particles is 
concerned : 

SU(4) : (T2/c7i = 1.3548 ±0.0064 adjoint : 1.3333 fundamental : 1.8182 
SU(6) : (T2/C71 = 1.6160 ± 0.0086 adjoint : 1.6000 fundamental : 1.9686 
(Ts/o-i = 1.808 ±0.025 adjoint : 1.8000 fundamental : 2.3635 

The results are that precise, that you see a two standard deviation from the adjoint, except 
for the second ratio of SU{6). This deviation is natural, since the diluteness of the magnetic 
quasi-particles is small, on the order of a couple of percent, as we will explain at the end of this 
subsection. So we expect corrections on that order to our ratios. 

'^Factorization is not obvious if uncorrclatcd Z(N) vortices'^ cause the area law: ~ {1 — cos{2kn/N)). This 
is due to the macroscopic size of the vortex perimeter causing long range correlations between loops. 
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There is a less preeise determination of the ratio 02/ (Ji = 1.52 ± 0.15 in SU{5) But the 
central value is within 1 to 2% of the predicted value 3/2 from the adjoint. The fundamental 
gives a ratio 1.8231. 

The SU{8) ratios are known on a rather course lattice ^'^ and using a different algorithm: 



c72/cri = 1.692(29) adjoint 
CT3/CT1 = 2.160(64) adjoint 
f74/c7i = 2.26(12) adjoint 



: 1.714 fundamental : 2.106 
: 2.143 fundamental : 2.958 
: 2.286 fundamental : 3.256 



In conclusion: the seven measured ratios are consistent with the quasi-particles being 
independent, as in a dilute gas and in the adjoint representation. The number of quasi-particle 
species contributing to the k-tcnsion is 2k{N — k). This number happens to coincide with the 
quadratic Casimir operator of the anti-symmetric representation. 

The fundamental monopoles are clearly disfavoured by the data. 

For other representations we have already said what to expect. E.g. the fully symmetric 
representation with k boxes will show the same tension, but through a very suppressed junction, 
especially for large N^. So lattice data will show predominantly the k{ai) tension. 

A caveat is in order. As we arc in three dimensions, the Coulomb law is logarithmic. It is 
multiplied by the quadratic Casimir C2(i?). So logarithmic confinement follows Casimir scaling! 

So if the fully symmetric representation with k boxes shows ratios compatible with the 
C2{Rsym) ^ k{N + k) the strings are probably behaving somewhere in between logarithmic 
and linear confinement. 

An inveterate pessimist could say the same of the antisymmetric reps where the Casimir 
ratio happens to be the same as for our adjoint monopole model. According to him, in the 
fitting of the tension the plateau has been mistaken for a linear law, whereas in reality it is 
logarithmic. I do not know whether the data exclude this. 

7 Conclusions 

There is remarkable agreement with an accuracy of about a percent between numerical sim- 
ulation at high T and the dilute adjoint monopole gas. This diluteness appears as a small 
parameter for every SU(N) theory. Some parameters, like the mass of the magnetic quasi- 
particle, are still within a large range, although its size is Im- Is it heavy, on the order of the 
lowest glueball mass (i.e. itlm), then the transition is that of a non-relativistic Bose gas. Is it 
light (with respect to its size) then the transition is that of a near-relativistic BE transition. 

Perhaps the best way to summarize our approach is to return to fig. 2. It is clear that the 
relation: 

a ~ lAiriM (28) 

implies that if^a is the diluteness 6, and must be small for consistency. And lattice data have 
borne that out! Also flux tube models like that of Isgur-Paton"^" predict this small ratio ^"'^ as a 
result of the balance of the string force and the phonon excitations of the closed string forming 
the lowest glueball. 

Now , once we accept the idea, that on the high temperature side of the transition a dilute 
gas describes things so well, the constancy of that diluteness down till T = suggests that all 

Tor SU(3) the k=2 symmetric representation still has an appreciable junction value, see ref. 
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what changes is the thermal wavelength Xb{T). At temperatures on the order of the ghieball 
mass or higher it is clear that the interparticle distance 0{l/g'^T) is larger than the thermal 
wave length 1/T, because the coupling is so small. But at a temperature Tq the thermal wave 
length takes over with Tq mmonopoie- It is very unlikely that Tq will be below Tc, since that 
would imply the BE transition as a second one below T^- So the monopole mass should not be 
lower than Tc, so the transition will then be a non- or near-relativistic BE transition. 

For the convenience of the reader we give the lattice numbers: 

Tc = (0.596(4) + 0.453(30)/Af2^cri and mQ++/ai = 3.341(76) + l.TS/iV^ 12 qj. j,^ = 
0.174(5)(1 + 0.27(3)/A''2)rno++. This gives the non-relativistic limit if the monopole is as heavy 
as rrimonopoie = Tno++ = Im- In the unlikely case that nimonopoie is as light as Tc we have the 
near-relativistic case. 

Below Tq the dilute gas gives a contribution to our ratios, but now determined by Bose 
statistics. For an analysis of the cold phase ratios see ref. 

Of course for most observables in the low T phase the Bose statistics is all-important. Thus 
the string tension will become non-zero below Tc, the character of the transition, i.e. a jump or 
continous behaviour in the occupation fraction of the p = states will be crucial to know and 
and is calculable in this model. Vortices in the condensed superfluid give then the flux tubes 
in the Isgur-Paton model of QCD. 

Realistic QCD involves quarks, and there is all reason to believe they couple strongly to 
our monopoles. After all, the latter are bound states of magnetic gluons. In the CFL phase 
colour is fully broken''^, and it is amusing to speculate on how the adjoint monopoles behave 
there. A fundamental multiplet of monopoles would get confined. 

It would be interesting to check this model in SUSY gluodynamics, where so many features 
are known analytically. 

The unbearable heaviness of the ground state (the energy density of our dilute gas is on 
the order of l{GeV)^) can be taken into account by string theory calculations of the same 
ratios Indeed, taking into account gravitational effects in an AdS/CFT context, does not 
affect the result gotten by the simple monopole gas picture: k{N — k) scaling does result in 
a large N limit, with k/N fixed. Note however that in practice these calculations are done in 
weak gravitational bulk fields. 

Note that a real time picture of our quasi-particles is lacking. According to our Euclidean 
picture they become, at very high T, a 3d gas of particles with small size Im and small inter- 
particle distance. The real time picture is a challenge. 

Luigi del Debbio, Pierre Giovannangeli, Christian Hoelbling, Dima Kharzeev, Alex Kovner, 
Mikko Laine, Biagio Lucini, Harvey Meyer, Rob Pisarski, and Mike Teper provided me with 
useful comments. 

I thank the organizers for their invitation, an inspiring meeting, and for wonderful hospi- 
tality. 
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